function [model, Ypred, blad, kor] = rownanieCsLProd(Y, U, na, nb, nc, nk)

%Y - LProd
%U - Capital Services

kor(1) = KorelacjaPrzes(U, Y, nk);
if (nb>1)
    kor(2) = KorelacjaPrzes(U,Y, nk+1);
end

DAT = iddata(Y, U, 1);
model = armax(DAT, 'na', na, 'nb', nb, 'nc', nc, 'nk', nk);

DATpred = predict(model, DAT, 1);
Ypred = DATpred.y;
blad = SumaKwadratowBl(Y, Ypred);